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ABSTRACT 

We present our new deep optical imaging and long-slit spectroscopy for Arp 220 that is the archetyp- 
ical ULIRG in the local universe. Our sensitive Ra imaging has newly revealed large-scale, Ka absorp- 
tion, i.e., post-starburst regions in this merger; one is found in the eastern superbubble and the other 
is in the two tidal tails that are clearly reveled in our deep optical imaging. The size of Ha absorption 
region in the eastern bubble is 5 kpc x 7.5 kpc and the observed Ha equivalent widths are ~ 2 A ± 
0.2 A. The sizes of the northern and southern Ha-absorption tidal tails are ~ 5 kpc x 10 kpc and ^ 
6 kpc x 20 kpc, respectively. The observed Ra equivalent widths range from 4 A to 7 A. In order 
to explain the presence of the two post-starburst tails, we suggest a possible multiple-merger scenario 
for Arp 220 in which two post-starburst disk-like structures merged into one, and then caused the two 
tails. This favors that Arp 220 is a multiple merging system composed of four or more galaxies, arising 
from a compact group of galaxies. Taking our new results into account, we discuss a star formation 
history in the last 1 Gyr in Arp 220. 

Subject headings: galaxies: individual (Arp 220) — galaxies: interactions — galaxies: starburst 



1. INTRODUCTION 

Ultra luminous infrared galaxies (ULIRGs) whose in- 
frared (8-1000 fim) luminosity exceeds 10 12 L© are the 
sites where extraordinarily active star formation occurs 
i n their central reg ions, triggered by galaxy mergers 
([Sanders et al.l [1988). Since such galaxies have been 
also found as massive, dusty galaxies at high redshift, 
providing substantial star formation in early universe 
(|Le Floc'h et al.l |2005[ ), it is important to understand 
their physical origins unambig u ously. 

Originally, ISanders et al.l (|1988[ ) proposed that 
ULIRGs come from a merger between two (gas-rich) 
galaxies. On one hand, multiple merger scenarios 
have also been proposed for the formation of ULIRGs 
based on their morpholog ical an d dynamical prop- 
erties (|T aniguchi_iLSIiio^ 119981 : iBorne et al.l 120001 : 
lHaan et al.ll201lf ). Unlike well-defined, on-going merging 
systems such as the Antennae (NGC 4038 + NGC 
4039), most of ULIRGs appear to have complicated 
morp hological properties of so -called advanced mergers 
(e.g., ISanders fc Mirabellll996D . 

Namely, the following questions have not yet been un- 
settled unambiguously for the origin of ULIRGs. (1) 
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How many galaxies were merged into one? (2) Which 
types of galaxies were merged into one? (3) How were 
their orbital parameters? Since the merging of galax- 
ies can smear out morphological properties of original 
galaxies that take part in the merger, it is generally dif- 
ficult to give firm answers to the above ques tions (e.g., 
IWeil fc Hernquistlll996l : iBournald et al.ll2QQ7D . 

It is worthwhile noting that a major merger be- 
tween two galaxies may be a site of multiple merg- 
ers because any gala xies have their satellite galaxies 
(|Wang fc Whi te 2012 and references therein). For ex- 
ample, let us imagine a case of the merger between 
our Galaxy and M 31. Since both galaxies have two 
significant satellite galaxies (LMC and SMC for our 
Galaxy and M32 and NGC 205 for M31), this major 
merger includes six galaxies. Since a minor merger 
also g ives dynamical disturbanc e to its host disk galaxy 
(e.g., [Mi nos & Hernquist 119961; [T aniguchi et lDXIl99ffl 
Bournald et al. 2007;"Villabos & Helmi 2008), it seems 
to have a unambiguous discrimination between a ma- 
jor merger between two galaxies and a multiple merger. 
It is also worthwhile noting that numerical simulations 
of multiple mergers can make an Arp 220-like morphol- 
ogy during the course o f mergers among six galaxies 
(Weil & Hernouiil [1996). A one-sided long tidal tail 
structure (e.g., Mrk 273) can als o be generated by a mul- 
tiple merger model (Bekki 2003). 

Among the ULIRGs in the local universe, Arp 220 is 
the archetypical advanced merger (jSanders et al.l 119881 : 
Uoseph fc Wri ght 1985|). Arp 220 shows two starburst 
events at least; one is the centra l super star clus- 
ters (i.e., the on-going starburst; iShaya et al.l Il994t 
iTaniguchi. Trentahm fc Shioval Il998h . and the other is 
the figure "8" sh aped structure of ionized gas d riven 
by a superwind ( iHeckman efall I1987L I1990L Il996f ). It 
is also suggested that there is a hidden starburst core 
in the central region of Arp 220 that is obscured seri- 
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ously , Ay ^30-50 mag (e.g., [A nantharam aiah et al.l 
2000; iShiova. Trentham. fc Taniguchil l2001f ). Indeed, 
multi-phase star formation events in Arp 220 and other 
ULIRGs have been identified in more recent imag- 
ing and spectroscopic investigations (e.g ., Wilson et al.l 
2006; [Rodriguez Zaurin et al.ll2QQ8L 120091 : ISoto fc Martini 



2010). Accordingly, Arp 220 appears to have experienced 



a long history of episodic starbursts during the course of 
merging within last ~ 1 Gyr. 

This complicated star formation history in Arp 220 
may come from a multiple mergers including more than 
two g alaxies (e.g. iTaniguchi fc Shioval[l998l : iBorne et al.l 
2000). It is generally difficult to explore the history 
of such multiple mergers in an advanced merger phase. 
However, since Arp 220 shows the evident on-going and 
past starburst events noted above together with a faint 
pair of tidal tails, it is worthwhile investigating its star 
formation history in more details. For this purpose, we 
have carried out deep Ha imaging using the Subaru Tele- 
scope and a long-slit spectroscopic observation with the 
Keck II Telescope. Based on our new observations, we 
discuss a possible scenario for the star formation history 
in Arp 220 and re- visit a multiple merger scenario. 

We use the distance of 77.6 Mpc toward Arp 220 that is 
estimated with both the syst emic velocity of 5434 km s _1 
(jde Vaucouleurs et al.|[l99lh and a Hubble constant Hq 
= 70 km s _1 Mpc -1 . At this distance, \" corresponds 
to 376 pc. 

2. OBSERVATIONS 
2.1. Optical Imaging 

In order to understand the whole star formation his- 
tory in Arp 220, one efficient method is to reveal post- 
starburst regions (their characteristic ages are ~ 0.2 - 1 
Gyr) as well as the on-going active star forming regions 
(their characteristic ages are < 10 Myr). For this pur- 
pose, we obtained deep Ha on-band and off-band frames 
using the Faint Object Camera And Spectrograph (FO- 
CAS) on the Subaru Telescope on 2002 September 15 
(UT). Total integration time for the on and off bands 
were 8400 seconds and 2400 seconds, respectively. The 
center wavelengths and bandwidths were 6709 A and 
85 A and 6588 A and 73 A for the on and off bands, 
respectively. The on-band frames include stellar contin- 
uum emission and Ha+[N II] AA6548,6583 line emission, 
while the off-band frames include stellar continuum only. 

Bias subtraction, flat-fielding, and sky subtraction 
were performed. Flux calibration was carrie d out by us- 
ing SDSS spectrum (jAbazajian et al.ll2009f ) of a galaxy 
in the same field of view. In continuum subtraction, we 
took consideration into the continuum slope between the 
Ha on and off wavelengths estimated from Keck spec- 
trum around eastern Ha absorption region (see Section 
[3]). This correction was applied to the whole region of the 
off-band image, and reduced continuum intensity by 2 % 
or Ha equivalent width by ~ 1.5 A than if we assume 
no continuum slope correction, or flat continuum. Since 
the eastern Ha absorption region and the two tidal tails 
have the same r' — i' color indices (Section [3]), this con- 
tinuum correction is reasonable. The spatial resolution 
of FOCAS images was 0".7. 

Deep i?-band frames of Arp 220 were also obtained 
by Subaru Prime Focus Camera (Suprime-Cam) on 2004 



Figure 1. (left panel) Continuum-subtracted, smoothed 
Ha+[N Il]AA6548,6583 intensity map, and (right panel) Ha equiv- 
alent width map of Arp 220 obtained by FOCAS. The color ranges 
are from —1 x 10 -18 to 5 x 10 -18 erg cm -2 s _1 (0.20 arcsec) -2 
for Ha+[N II] intensity map, and from 10 to —30 A (positive value 
means absorption) for Ha equivalent width map. The contours are 
shown at 2, 6, and 10 A in the Ha equivalent width in absorption. 
North is up and east is left. The slit position in the longslit spectro- 
scopic observation is shown with the diagonal green box. Two red 
rectangles represent the positions where continuum surface bright- 
ness and equivalent widths are shown in Figure [3] Green bar in 
each panel at lower right corresponds to 10 kpc. 

February 18 (UT). Bias subtraction, flat-fielding, and sky 
subtraction were performed. Total integration time was 
4200 seconds, and the spatial resolution was 1".0. This 
i?-band image was used to find faint tidal features around 
Arp 220 (see also. lKoda et al.ll2009[ ). 

2.2. Optical Spectroscopy 

A new optical spectrum was obtained by the Low Res- 
olution Imaging Spectrometer (LRIS: lOke et afl [l995) 
with the Keck II Telescope on 1997 April 16 (UT). A 
900-second exposure was taken with the 600 line mm -1 
grating, with a central wavelength setting of 6024 A. The 
grating is blazed at 5000 A and the dispersion is 1.28 A 
pixel -1 . The observed wavelength range and the spectral 
resolution were 4670 - 7210 A and 6.9 A, respectively. 
The 1" slit with a length of 66" was oriented at the 
position angle (PA) of 130 degree to observe both the 
nucleus of Arp 220 and its visual companion D (seen at 
SE of Arp 220) simultaneously. The spectrum was re- 
duced by using IRAF with the standard manner, i.e., 
wav elength calibration, f lat-fielding, and flux calibration 
(see lOhvama et aflll999l for detail). 

3. RESULTS 
3.1. Imaging Data 

In post starburst regions, massive OB stars already 
ceased and thus intermediate-mass stars such as A- 
type stars are dominant populations. Such regions are 
probed by strong hydr ogen Balmer absorption lines (e.g., 
ITaniguchi et all 119961 ). 

Figure [1] shows the continuum-subtracted, smoothed 
Ha+[N II]AA6548,6583 image, and the Ha equivalent 
width maps obtained by FOCAS, respectively. We detect 
strong Ha emission from the central region of Arp 220 
and shell-like structures in both SE and NW directions, 
as already repo rted in the literature ([Heckman et al.l 
1987. Il990l [1996) . The strong Ha emission from the cen- 
tral region is considered to be attributed to the on-going 
starburst activity there. On the other hand, the shell-like 
structures is due to the recent-past superwind, which is 
a galactic-scale outflow caused by many supernovae. 

On closer inspection, we find three Ha absorption re- 
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Figure 2. i^-band image of Arp 220 obtained by 
Suprime-Cam (right panel) with the post starburst re- 
gions shown by red contours. For comparison, the 
HST/ACS ima ge of Arp 220 is also show n in t he left panel 
(http://hubblesite.org/gallery/album/pr2008016aq/ ). North is up 
and east is left. The post starburst regions are shown by red-color 
contours. Green bar corresponds to 10 kpc. 
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Figure 3. Continuum surface brightness and Ha equivalent width 
profiles at (a) northern and (b) southern tails using FOCAS data. 
Solid lines represent continuum surface brightness in arbitrary unit 
(—10 A corresponds to no flux), while filled circles represent equiv- 
alent widths. The silt positions are shown in Figure^ The aper- 
ture size in each data point is 2" x 8 /x . Uncertainties in continuum 
surface brightness are almost the same at all positions, and only 
the rightmost one is shown. 

gions. One Ha absorption region is detected about 20" 
east from the center, inside the eastern Ha bubble. The 
observed equivalent width of Ha in the center 3" x 3" 
is 2.3 A ± 0.15 A (la). The size of this region is ap- 
proximately 5 kpc x 7.5 kpc. The other two spectacular 
Ha absorption regions appear to be associated with the 
two optically-faint tidal tails found in our deep i?-band 
image (see also Figure [2]). The observed equivalent width 
of Ha amounts to ~ 5 ± 0.4 A in total; ~ 4 ± 1.1 A at 
the continuum flux peak, and ~ 7 ± 1.7 A at the outer 
side of the tails (Figure [3]). The sizes of the northern and 
southern Ha absorption regions are ~ 5 kpc x 10 kpc 
and ~ 6 kpc x 20 kpc, respectively. Since these three 
regions show the Ha absorption feature, A- type stars, 
or post-starburst populations, are dominant in the opti- 
cal continuum without Ha emission sources such as OB 
stars. 




Figure 4. (a) u' — g f , (b) g' — r' , (c) r' — i' , and (d) i' — z' color 
index maps obtained from SDSS images. The color scales are from 
— 1 to 3 mag for u' — g' map, while from —1 to 1 mag for the other 
maps. North is up and east is left. White bar in each panel at 
lower right corresponds to 10 kpc. 

Figure [4] displays Arp 220 color index maps obtained 
from SDSS images (jAbazaiian et al.ll2009D . We find that 
the three Ha absorption regions have almost the same 
color indices. Furthermore, their color indices within 
several kpc except for the nuclear region are similar to 
those of the Ha absorption regions. These facts suggest 
that the dominant stellar population at these regions 
are t he same, i.e., post-star burst, which is consistent 
with Rodrig uez Zaurm et al.l (j2008[ ) and iSoto fc Martini 
(2010) results. 

Finally, we comment on an alternative interpretation 
that Ha absorption in the two tidal tails could be caused 
by t he high extinction found in a number of ULIRGs 
(e.g., Rodrig uez Zaurm et al.l [20091 ). If the extinction is 
very high enough to hide on-going massive star forma- 
tion, the observed optical light might be dominated by 
intermediate- type stars, resulting in strong Ha absorp- 
tion with little emission line. Although such situations 
may occur in the central parts of ULIRGs including Arp 
220 jSh lova, Trentham. fc Taniguc hi 2001), it seems un- 
likely that such very dusty starbursts occurred in the en- 
tire regions of both the eastern bubble and the two tidal 
tails. Therefore, we consider that the post-starburst in- 
terpretation is a more robust idea for the observed Ha 
absorption regions in Arp 220. 

3.2. Spectroscopic Data 

The so-called Lick index (jWorthev fc Ottavianil Il997l ) 
provides a common method for calculating an equiva- 
lent width. We therefore calculated the H/3 Lick index 
as a function of slit position. Ha equivalent width was 
also calculated in the same way as that of Lick index. 
Note that we adopt the following wavelength coverage; 
6552.4 - 6564.6 A for the on-Ha region, 6505.9 - 6525.44 
A and 6603.7 - 6664.9 A for the off-Ha region. We 
chose these ranges in order to avoid the effects from 
[N II]AA6548,6583 emission from the object, other ab- 
sorption lines, and sky emission. 

Figure [5] represents continuum flux density and equiva- 
lent width profiles along the slit position displayed in Fig- 
ure [TJ The Ha equivalent width at the southern part of 
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Figure 5. Continuum flux density and equivalent width profiles 
along the slit displayed in Figure [TJ using LRIS spectroscopic data. 
The horizontal axis represents the distance from the continuum 
peak in optical wavelength. Positive value in horizontal axis means 
in the northwest of the nucleus. Solid line represents the surface 
brightness of continuum in arbitrary units (—10 A corresponds to 
no flux). Filled circles and open triangles represent the Ha equiv- 
alent widths and H/3 Lick index, respectively. Positive equivalent 
widths mean absorption, while negative ones emission. 

the eastern bubble is ~ 1.5 A, i.e., Ha absorption, which 
is consistent with FOCAS equivalent width map (Figure 
[3]). The H/3 equivalent width is more sensitive to post- 
starburst population than that of Ha, because the flux 
of the H/3 emission lin e is generally about 3 times s maller 
than the Ha flux (e.g. JOsterbrock fc Ferlandll2QQ6[ ). The 
observed Hf3 Lick index of the eastern Ha absorption 
region is 3 - 4 A. H/3 absorption is detected not only 
around eastern bubble but also in the central region of 
Arp 220, ~ +5" and ~ — 5" from the center (Figure 
EJ. This fact suggests that post-starburst populations 
spread to several kpc in diameter, which is consistent 
with the results of b oth [Rodriguez Zaurm et al.l {2008) 
and iSoto fc Martini (|2010l ). Unfortunately, we cannot 
measure the H/3 equivalent width of the west Ha ab- 
sorption tails, because they were out of the field of view 
of LRIS. 

4. DISCUSSION 
4.1. A Multiple Merger Scenario 

Based on its morphology, it is obvious that Arp 
220 is now in an a dvanced phase of a galaxy merger 
([Sanders et all 119881 ). For Arp 220. ITaniguchi fc Shioval 
(| 19981 ) have discussed a possibility of multiple mergers. 
Their argument is based on these facts that Arp 220 
has two compact starburst regions and that a pair of 
nuclei (i.e., a pair of supermassive black holes) is nec- 
essary to initiate the so-called nuclear starburst in a 
merging system in eithe r a minor or a major merger 
(Taniguchi & Wada 1996). This may be a plausible idea 
because the starburst is intrinsically different from ord i- 
nary star formation in galactic disks (Dad di et al.l l2010). 
The presence of multiple (at least three) OH megamaser 
sources can be interpreted as evi dence for major merg - 
ers of four galaxies in Arp 220 (Dia mond et al.| [l989). 
However, although further observations detected multi- 
ple components, it is reminded t hat they do not nec es- 
sarily suggest the multiple nuclei ([Rovilos et al.l 12003ft . 

One int eresting advantage of the multiple merger 
model by ITaniguchi fc Shioval (|1998[ ) is that they pre- 
dicted the presence of a pair of counter rotating dense 
molecular gas disks before t he ac tual observational de- 
tection by 1 Sakamoto et al.l (|1999). Arp 220 is a very 



gas rich system with ~ 1O 1O M and a half of gas is as- 
sociated with two a pparent eastern and western nuclei 
(|Scoville et al.l 1199/1 ). It seems difficult to estimate or- 
bital parameters of merging galaxies to explain all these 
observational results. 

However, the post-starburst tidal tails revealed by our 
our new observations presented here suggest a possible 
new idea on the merging history in Arp 220. For further 
consideration of multiple mergers, we discuss the origin 
of Arp 220. 

Here we focus on the origin of two post-starburst 
tidal tails. First, in antenna-like tidal tails, one tidal 
tail emanates from one colliding gal axy and thus two 
well- d eveloped tidal tails are formed (jToomre fc Toomrel 
Il972[ ). Taking this into account, we can postulate that 
two post-starburst systems merged into one and then 
the observed two post-starburst tails were formed in Arp 
220. The observed smooth distribution of Ha absorp- 
tion features requires this scenario although small-scale 
star formation activity could occur in gaseous tid al tails 
(jBarnes fc Hernquistlfl99l iDuc fc Mirabellll998[ ). 

It should be noted that the size of each post-starburst 
system must be as large as several kpc at least, corre- 
sponding to a typical scale of disk galaxies. Let us con- 
sider a probable star formation history in a major merger 
between two disk galaxies. During the course of such 
merger evolution, enhanced star formation events oc- 
curred, being triggered by inward gas flo w driven by non- 
axisymmetries in the galaxy disks (e.g., Di Ma tteo et al.l 
2007). However, it should be reminded that the star- 
burst events are essentially confined in the nuclear re- 
gion of each galaxy; i.e., within only a few kpc diameter 
area. Note that this is comparable to t hose of nuclea r 
starburst regions in disk galaxies (e.g., Balzano 1983). 
In order to make a pair of very large-scale (~ 10 kpc), 
post-starburst tidal tails observed in Arp 220, we need 
two large-scale (~ several kpc) post-starburst disk-like 
structures. Such a structure can be made in an advanced 
phase of a merger between two disk galaxies; i.e., a ma- 
jor merger remnant. We therefore strongly suggest that 
the observed two tidal tails in Arp 220 need a merger 
between two advanced (i.e., post-starburst) merger rem- 
nants. Namely, we need four disk galaxies to explain the 
observed post-starburst tidal tails in Arp 220. We then 
conclude that Arp 220 comes not from a typical merger 
between two galaxies but from a multiple merger includ- 
ing four galaxies at least. 

We estimate the travelling time of stars in the tidal 
tails from the merger center to the present place of tails. 
Given both the velocity of tidally-liberated stars as 100 
km s _1 and the distance as 20 kpc, the travelling time is 
calculated as 200 Myr, being consistent with the lifetime 
of post-starburst phase (0.2 - 1 Gyr). Our scenario for 
the evolution of Arp 220 is summarized in Figure [6] 

4.2. Star Formation History in Arp 220 

Finally, we consider a possible unified picture for var- 
ious star formation events in Arp 220. As shown in 
Figure 1, there are two more important starburst fea- 
tures in Arp 220; one is the central starburst and the 
other is t he pair of superbubble caused by the superwind 
activ ity (jHeckman et al.l 119871 11990L 119961 : IColina et all 
2004). The central starburst must be occurred very 
recently, < 10 Myr ago, being the life time of mas- 
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Figure 6. A possible scenario for the formation of the post- 
starburst tails. (1) Two merging partners caused an intense star- 
burst. Here we postulate that a couple of intense starbursts (SB) 
occurred in the merging system about 0.2-1 Gyr ago. (2) After the 
cease of the starburst, these regions evolved into the post-starburst 
(post-SB) regions. And, (3) they merged into one and then the two 
post-SB tails were formed. 




Figure 7. A possible unified picture for various star formation 
activities in Arp 220. 

sive OB stars taken into account. On the other hand, 
since the dynamical timescale of the obs erved superbub- 
ble is estimated as a few tens Myr (e.g., iHeckman et al.l 
1990), Arp 220 experienced another starburst a few 
tens Myr ago. These star formation events corres pond 
to those found in IRodrfguez Zaurm et al.l ((2008) and 
[Soto & Martin (2 0ToF i.e., the young stellar population 
in IRodrfguez Zaurm et al.l (|2008[ ) and the const ant- star- 
formatio n region occurring in the central kpc region of 
Arp 220 (ISoto fc Martini l201(l . 

In summary, taking account of the two post-starburst 
tails, we propose a unified evolution model for Arp 220 
as shown in Figure [3 These three starburst events could 
be driven by sequential multiple mergers in Arp 220. 

5. CONCLUDING REMARKS 

There are a significant nu mber of compact groups of 
galaxies in the local universe (Hickson 1982ab), leading 
to multiple mergers. Such compact groups of galaxies are 
also present in rich groups and out s kirt regions of cluste r 
of galaxies at intermediate redshift (|Mendel et al.l l2011). 
Moreover, a possible multiple merging system has been 



also found at z > 1 ([Schawinski et al.ll2011l ). 

ULIRGs are considered to evolve to quasars and then 
giant early type galaxies (elliptical or SO galaxies). A 
large number of high-redshift analogs of ULIRGs, so- 
called submm gal axies (SMGs) have been discovered be- 
yond z ~ 1 fe.g.JHughes et al.l 1998; Bargerel_al][l998; 
lAmblard et alJl201lUMagdis et al.ll201lD . Although ma- 



jor mergers between only two galaxies have been often 
considered in terms of the hierarchical evolution of galax- 
ies from high redshift to the present day, we will have to 
take the impact of multiple mergers into account in fu- 
ture consideration. 
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